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Abstract 

Background: Tuberculosis remains a high burden for Human society despite considerable investments in its 
control. Unique features in the history of infection and transmission dynamics of tuberculosis pose serious 
limitations on the direct interpretation of surveillance data and call for models that incorporate latent processes and 
simulate specific interventions. 

Methods: A transmission model was adjusted to the dataset of active tuberculosis cases reported in Portugal 
between 2002 and 2009. We estimated key transmission parameters from the data (i.e. time to diagnosis, treatment 
length, default proportion, proportion of pulmonary TB cases). Using the adjusted model to the Portuguese case, 
we estimated the total burden of tuberculosis in Portugal. We further performed sensitivity analysis to 
heterogeneities in susceptibility to infection and exposure intensity. 

Results: We calculated a mean time to diagnose of 2.81 months and treatment length of 8.80 months in Portugal. 
The proportion defaulting treatment was calculated as 0.04 and the proportion of pulmonary cases as 0.75. Using 
these values, we estimated a TB burden of 1.6 million infected persons, corresponding to more than 15% of the 
Portuguese population. We further described the sensitivity of these estimates to heterogeneity. 

Conclusions: We showed that the model reproduces well the observed dynamics of the Portuguese data, thus 
demonstrating its adequacy for devising control strategies for TB and predicting the effects of interventions. 
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Background 

Twenty years after considered to be a major global pub- 
lic health problem, tuberculosis (TB) remains a high 
burden for human society. In fact, WHO places TB as 
the second leading cause of death from an infectious 
disease, with almost two million deaths per year [1]. 
Moreover, according to this organization, one third of 
the world's population is infected with TB and nearly 10 
million people progress to an active state of the disease 
every year. Some aspects of TB transmission are slower 
than others in responding to changes in biologic, demo- 
graphic, and socioeconomic conditions, resulting in a 
superposition of time scales that poses major challenges 



* Correspondence: j. sollari. lopesiagmail.com 

'institute Gulbenkian de Ciencia, Apartado 14, 2781-901 Oeiras, Portugal 
Full list of author information is available at the end of the article 

Bio Med Central 



to the interpretation of transmission measures. A frac- 
tion of infected people progresses to active disease 
within two years of infection, while others may maintain 
a latent infection for decades with a reduced risk of dis- 
ease progression. Complicating things further, those who 
have had the disease and have been successfully cured also 
maintain a risk of recurrent disease [2]. Furthermore, host 
susceptibility and rates of activation and reactivation of 
the disease appear to be related to co-morbidities (e.g. 
malnutrition, alcoholism, drug-addiction, diabetes, or 
HIV infection and other immunosuppressive diseases), 
high-risk of exposure to infection (e.g. health workers, 
prisoners), ageing and socioeconomic factors [3,4]. A 
challenge for mathematical modeling is to capture those 
aspects that are essential to the questions being ad- 
dressed while maintaining a stylized structure that enables 
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interpretation without compromising accuracy of the 
conclusions [5-8]. In this respect, the selection of the 
model is intimately connected to the posed question [7], 
but also to specificities of the settings being studied. 

In this work we applied a mathematical model of TB 
transmission to Portuguese surveillance data. Recorded 
data on TB incidence in Portugal from ECDC [9] show a 
consistent decrease from 1999 to 2011 of 4.34% per year, 
resulting in a total decline in TB incidence of 52.13%. 
Nevertheless, Portugal registers the highest incidence 
among Western European countries (over 20/100,000, 
while a great majority is below 10/100,000), which re- 
flects a considerable rate of TB transmission. Recent 
immigrants comprise only about 3% of overall yearly 
cases, and the 50% most common ages of TB patients 
vary vastly from 25 to 65 years old (note also that TB 
cases in children with less than 15 years old comprises 
an uncommonly significant proportion of almost 3%) 
[10]. Other common risk factors have some weight on 
the overall yearly TB cases in Portugal, for example being 
a prison inmate (about 2%), HIV-positive (about 14%) 
or drug-user (about 15%) [10]. Although not modeled 
explicitly, these factors are generally contemplated by 
our analysis of sensitivity to risk heterogeneity. 

Initially, we applied a mathematical model for TB trans- 
mission under previously specified homogeneity assump- 
tions [11]. For the analysis, we used data records from 
2002 to 2009 provided by National Directorate of Health 
(Additional file 1). Using these data, we estimated param- 
eter values of the model (i.e., detection time, treatment 
length, proportion of unsuccessful treatment, and propor- 
tion of pulmonary TB cases). Subsequently, we assessed 
the fit of the model by contrasting dynamic trajectories de- 
scribed by the parameterized model with real dynamics 
described by the data. The model was used to calculate the 
total burden of TB in Portugal considering the different 
stages of infection (i.e., primary infection, active infection 
not yet diagnosed, latent infection, and active infection 
under treatment). We further consider the existence of dif- 
ferential susceptibility to infection by applying an extended 
heterogeneous model based on previous work [8] . 

Methods 

Analysis of treatment length and outcome in Portugal 
from 2002 to 2009 

Surveillance data on TB for the period of 2002-2009 was 
provided by the Portuguese National Directorate of Health. 
These data, collected by medical practitioners, consisted of 
date of first symptoms, date of diagnosis and treatment 
starting, date of treatment ending, clinical form (i.e. pul- 
monary or extra-pulmonary), and treatment outcome [i.e. 
treatment completed (with or without laboratory confirm- 
ation); treatment failure; treatment default (i.e. lost to fol- 
low up); death (from tuberculosis or other causes)]. These 



data had been previously anonymized and, thus, no ethical 
approval was needed for the analysis presented here. From 
the dates of first symptoms and treatment starting and 
ending, we calculated time to detection and treatment 
length for each reported TB case - while the latter calcula- 
tion is clearly objective, the former can raise some ques- 
tions since retrospective self-report of symptoms is only a 
rough approximation for infection time. The proportion 
of treatment incompleteness was calculated by pooling 
together cases of default and treatment failure. The time 
to detection (1/r yrs), treatment length (1/oVyrs), pro- 
portion of incomplete treatments {(p T ) and proportion 
of pulmonary TB cases (v) were estimated directly from 
the Portuguese dataset. To accommodate possible outliers 
from reporting/clerical errors or atypical cases, we dis- 
carded values of treatment length and detection time lon- 
ger than 3 yrs. Following parameter estimations from raw 
data, cumulative data sets were created. Goodness-of-fit be- 
tween the cumulative data and the theoretical expectation 
was assessed using two-sample Kolmogorov-Smirnov tests. 

Basic TB transmission model 

In order to analyze the dynamics of TB epidemiology, we 
considered a TB transmission model modified from 
Gomes et al. [11]. In this model, susceptible individuals 
(class S) become primary infected with TB (class P) at rate 
A yrs 1 . From this class, a proportion <p enters an actively 
infected (and infectious) state (class I) at rate r5 yrs~ , 
whereas the remaining individuals enter a latent infected 
state (class L) at the same rate. Actively infected individ- 
uals start TB treatment (class T) at rate tk yrs -1 , where r 
is the inverse of the average time to detection and k is the 
proportion of active TB cases that is actually detected and 
initiates treatment per year. Individuals under treatment 
are assumed to be neither infectious nor susceptible to re- 
infection. These individuals leave treatment at a rate S T 
yrs 1 . A proportion (p T leaves class T to I due to either 
treatment failure or default, while the remaining (1 - <pr) 
are successfully treated and transferred to class L. The 
model is described by the following system of ordinary dif- 
ferential equations: 

dS . . 

— = fi-(X + fi)S 

dP 

— = AS + o\L-{8 + u)P 
at 

— =(f>5P+0)L + (f> T S T T-(TK + l4.)I (1) 

at 

= (l-4»)8P + {l-4> T )S T T-{aX + co + fi)L 

dr 

— = tkI-(S t + fi)T 
at 

The model parameters along with their typical values 
used herein are listed in Table 1. The force of infection 
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Table 1 Parameters of tuberculosis transmission models 



Symbol 


Definition 


Value 


(3 


Transmission coefficient 


variable (yrs" 1 ) 


V 


Proportion of pulmonary TB cases 


0.75 


I" 


Death and Birth rate 


1/80 yrs" 1 


6 


Rate at which individuals leave 
P compartment 


2 yrs 


<P 


Fraction of infected population 
developing active TB 


U.Ud 


a 


Reinfection factor 


0.5 


(i) 


Rate of endogenous reactivation 


0.0003 yrs" 1 


T 


Inverse of time to detection 


4.26 yrs" 1 


K 


Proportion of detected cases in a year 


0.87 


6 T 


Inverse of treatment length 


1.36 yrs" 1 


<Pt 


Fraction of treatment default and failure 


0.04 


(f 


Low-risk factor 


0.15 or variable 


f 


Proportion of low-risk group 


0.98 or variable 



a - Parameter used only in the heterogeneous model. 



is simply X = [ivl, where /> is the transmission coefficient 
and v is the proportion of pulmonary TB cases. Birth and 
death rates are assumed equal, fixed as \i = 1/80 yrs -1 . The 
reinfection factor a is the decreased chance of contriving 
TB by exogenous reinfection after previous exposure (ei- 
ther after active disease followed by successful treatment 
or after having only remained in a latent state). For this 
factor we considered a= 0.5, as estimated in Gomes et al. 
[8] using a world-wide dataset. For the endogenous reacti- 
vation rate co, we adopted cj = 0.0003 yrs 1 , as estimated 
by Vynnycky et al. [12,13] for the population of the United 
Kingdom. 

Parameters describing the transition between the pri- 
mary infection class P to L and / were revised. Numer- 
ous TB models do not consider a primary infection 
class P explicitly (see [7] and references therein). In our 
model, this simplification would correspond to setting 
the rate of progression out of P as 8 — > °°. Instead, we 
specified a primary infection stage based on information 
from the literature. After invasion of TB bacteria, some 
balance between the host immune system and the 
pathogen occurs. In about 5% of cases, the infection 
progresses to an active state within 2 years, and an add- 
itional 5% become actively infected later on [2,14]. 
Hence, we specified a rate (S) out of the primary stage, 
P, and a proportion (</>) to represent those progressing 
directly from primary infection to active disease. As- 
suming that 5% of the individuals become actively in- 
fected within 2 years of exposure, we set parameter 
values compatible with observed rates, as <p = 0-05 and 
5=2 yrs" 1 (see Figure 1). Treatment class is described 
by four parameters, specifically rate of detection (r), pro- 
portion of detected cases (k), treatment length (l/d T ) and 




0.0 0.5 1.0 1.5 2.0 2.5 

time (yrs) 



Figure 1 Sojourn time distribution in the primary infection 

state. Exponential distribution for the time infected individuals 

remain in primary state (P), and proportional progression to active 

disease (I) and latent (L) states. The curves are given by P(f) = e" a ' ', 

l(t) = 4>- (1 -e~ 6t ), 1(f) = (1 -<j>)- (1 -e" 6 -'), with parameter values 

0 = 0.05 and 6 = 2 yrs" 1 , in agreement with the expectation that 5% 

of primary cases progress to active disease within 2 years of 

infection (dotted-line). 
^ ) 

the proportion of incomplete treatment (</> r ). Parameter k 
was taken from WHO estimates for Portugal as 87% [1]. 
The remaining parameters were estimated directly from 
the Portuguese dataset as stated before. 

The basic reproduction number R 0 was obtained using 
the approach in [15]: 

R = PvSVt + pMh + o) 

fi(8 + fi)[(ji + cj)(tk+ S t + ft) + S t tk(1-4> t )] 

(2) 

System (1) has two equilibria: a disease-free equilib- 
rium Ei = (1,0,0,0,0) stable when R 0 < 1; and an endemic 
equilibrium E 2 = {S ,P ,1 ,L ,T ), where all coordinates are 
nonzero, stable when R 0 > 1. Details on the model and 
its equilibrium points are in Additional file 2. 

TB model assuming heterogeneous infection risk 

The mathematical model represented by system (1) can 
be extended as in [8] to enable a heterogeneous risk of 
infection in the population. In this implementation the 
population is composed of a subpopulation 1 with low 
risk of infection, where the force of infection is affected 
by a factor a.\ < 1, and a subpopulation 2 with high risk 
of infection, where the force of infection is increased by 
a factor ai > 1. Parameters yi an d Yi are t ne proportions 
of the population in each group, such that y 1 + y 2 = 1. 
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The average risk factor is normalized by considering 
Y\a.Y + X2«2 = 1» thus, heterogeneity is fully parameterized 
by cti {=a) and y 1 (=j). The model is expressed as 

-~ = Yffi-(Xi + fi)Si 

dP 

-± = \ i S i + o\ i L i -{8 + i*)P i 
at 

= (frdPi + 0)Li + (f> T S T Ti-(TK + fi)Ii 

at 

^ = (l-4>)8P t + (l-^SrT-iaXi + u + ft)L t 

= TKli-(S T + ft)Ti 

(3) 

where A,- = afivl, i = {1,2} and a and y can take values from 
the interval (0,1). For reference, [8] estimated a = 0.15 and 
y = 0.98 using a world-wide dataset. 

The basic reproduction number R 0 for system (3) is also 
expressed by equation (2). This can be derived by noting 
that in the heterogeneous model the average number of 
infections contributed by a single infectious individual 
is the sum R 0 = y^iRo + y 2 <X2Ro , where the first and 
second terms account for transmission to the low- and 
high-risk groups, respectively. Considering that y\ + y 2 = 1 
and y 1 a l + y 2 a 2 = 1, and that y\ = y and a.\ = a, we obtain 
R Q = yaRo + (l-ay)Ro, or simply R Q = Rq. From now on 
the basic reproduction number is represented by R 0 ir- 
respective of whether we refer to the homogeneous or 
heterogeneous model. 

System (3) has also two equilibria, a disease-free equi- 
librium E[ = (y,l - ^,0,0,0,0,0,0,0,0) stable when 2? 0 < 1, 
and an endemic equilibrium E 2 = (SJ, S 2 , P\, P* 2 , I\, I 2 , 
L\, L* 2 , T*, T 2 ), whose coordinates are all nonzero for 
the considered range of parameter values, which is stable 
when R 0 > 1. 

Results 

Parameter estimation from TB data in Portugal 

For the analysis on detection time, we selected 16,257 
cases (where detection time < 3 yrs) corresponding to 
99.7% of all data. The calculated mean value of detec- 
tion time and its standard error was 2.81 ± 0.02 months 
(r = 4.26 yrs -1 ), while its median was 2 with an interquartile- 
range also of 2. Given these parameters, the model ex- 
pectation for the flow into the under-treatment class T 
of a cohort I 0 is 

/o-r(t) :=I 0 {l-e- TK t ) (4) 

where I 0 = 16,257, r = 4.26 yrs 1 and k = 0.87. Figure 2a 
shows a superposition of expression (4) and the cumula- 
tive frequency of detection times of the selected cases. 
The agreement between the exponential distribution and 



the observed sojourn times in the infectious class / is 
impressive given that only the central tendency of the 
raw data was fitted (D = 0.243, p-value = 0.226). 

For the analysis of treatment length we selected 
17,477 cases (where treatment length < 3 yrs) corre- 
sponding to 99.9% of all data. The calculated mean 
value of treatment length and its standard error was 
8.80 ± 0.03 months (5r=l-36 yrs 1 ), while its median 
was 9 with an interquartile-range of 4. As before, the 
model expectation for the flow out of the under-treatment 
class T of a cohort T 0 becomes 

To^I(t)+L(t) := r 0 .(l- e -^" £ ) (5) 

where To = 17,477. Figure 2b shows a superposition of 
expression (5) with S T = 1.36 yrs 1 , and the cumulative 
frequency of treatment length of the selected cases. These 
results also show some agreement between the observa- 
tions and the exponential curve (D = 0.4595, p-value = 
6.746 x 10~ 4 ). 

As for the ratio of treatment incompleteness (i.e. treat- 
ment default or failure), we considered the data from 
2002-2009 with known outcome. The ratio of defaults 
and failures to the total number of treatments was calcu- 
lated as 0.04. Note also that, in the case of the Portuguese 
dataset, the rate of failures is insignificant compared to the 
rate of defaults (Table 2), thus, <p T can be approximated as 
the rate of defaults in TB treatment outcomes. Consider- 
ing the flow out of T 0 to class / only, the theoretical recon- 
struction can be written as 

IW(f) := T 0 -{l-e- ST t )-<p T (6) 

where T 0 = 17,477. Figure 2c shows the superposition of 
expression (6), where St= 1-36 yrs~ and (pr= 0.04, and 
the cumulative frequency of treatment incompleteness 
for the selected cases. Again, there is an impressive 
agreement between the observations and the theoretical 
curve, showing that sojourn times in treatment Tare also 
approximately exponential (D = 0.3226, p-value = 0.079. 
Ties on observed data were discarded for a correct calcu- 
lation of the p-value). 

Finally, for the proportion of pulmonary TB cases, we 
considered once again the full data from 2002-2009 
(Table 3). The ratio of pulmonary TB cases was 0.75. 

Epidemiology of TB in Portugal 

ECDC recorded data show a consistent decline from 
1999 to 2011 on TB incidence in Portugal [9], suggesting 
a non-steady state of TB transmission. Nevertheless, as 
an initial approximation, this study assumes equilibrium 
conditions for TB transmission, with an average number 
of reported TB cases of 2348.25 cases per year (Table 2). 
According to the Portuguese National Institute of Statis- 
tics, the population size of Portugal during this period 
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r \ 

Figure 2 Superposition of observed data and theoretical 
expectations, a) Cumulative frequency of detection of tuberculosis 
infected patients and theoretical expectations for the inflow to class 
fassuming t = 4.26 yrs , b) Cumulative frequency of TB treatment 
length and theoretical expectations for the outflow from class T 
assuming 6 T = 1.36 yrs -1 . c) Cumulative frequency of TB treatment 
defaulters and theoretical expectations for the flow from class T to 
class / assuming 67-= 1.36 yrs - ' and 0 7 =O.O4. 



was approximately 10.56 million people, resulting in a 
rate of 22.24 treatments initiated per 100,000 person- 
years. 

From system (1) and considering that TB transmission 
is at equilibrium, TB incidence can be calculated as 

Y = <pSP + oL + <p T S T T = {tk + p)I (7) 

considering that tk > > fi, we can simplify expression (7) as 
Y= tkI. Thus, TB incidence in Portugal can be simply cal- 
culated as the rate of entering treatment Y= 2.224 x 10~ 
per person-years (i.e. 22.24 per 100,000 person-years), 
while the prevalence is approximated by / = YItk = 
4.261 x 10~ 5 (i.e. 4.261 per 100,000 persons). Assuming 
that system (1) is at equilibrium and considering the 
parameter values of Table 1 and that I = 4.261 x 10~ , 
we calculated the value of ft and of the proportion of in- 
dividuals in the remaining classes of the model (Table 4). 
Thus, we obtained R 0 = 1.116 (yS = 72.358 yrs" 1 ). These 
estimates resulted in expectations for T that are congru- 
ent with the data collected in Portugal (see Additional 
file 2). Sensitivity analysis to u, a and S showed the ro- 
bustness of these estimates (Additional file 3). Note that 
given the uncertainties regarding TB reinfection, we con- 
sidered a range of values for sensitivity analyses for a that 
encompass the possibility that a previous TB infection 
increases risk of getting reinfected [16]. 

Considering the Portuguese population size, the esti- 
mated values correspond approximately to 11,000 people 
with primary TB infection, 1,600,000 with latent infection, 
1,200 patients under treatment and 500 undiagnosed 
actively infected persons. 



Table 2 Outcome of tuberculosis treatment in Portugal, 
2002-2009 



Outcome 


2002 


2003 


2004 


2005 


2006 


2007 


2008 


200« 


Death 


126 


126 


145 


142 


129 


130 


122 


100 


Default 


110 


93 


92 


98 


102 


61 


65 


46 


Failure 


1 


6 


3 


3 


4 


4 


0 


0 


Completed 


2,109 


2,074 


2,162 


2,091 


2,146 


2,045 


1,971 


1,1 9i 


Unassigned 


44 


39 


50 


26 


42 


91 


125 


865 
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Table 3 Clinical form under tuberculosis treatment in Portugal, 2002-2009 

Clinical form 2002 2003 2004 2005 2006 2007 2008 2009 

Pulmonary 1,748 1,717 1,844 1,752 1,850 1,749 1,745 1,643 

Extra-pulmonary 643 621 610 608 573 582 538 566 



Effect of considering a heterogeneous model 

Following the same reasoning as used for system (1), we 
assumed for system (3) that TB transmission in Portugal 
from 2002 to 2009 is close to equilibrium and that the 
fraction of actively TB infected individuals (/J + /J) was 
4.261 x 10~ 5 . By considering the parameter values of 
Table 1, we obtained estimates of the stationary values 
(Table 4) and R 0 = 2.237 (yS = 145.024 yrs" 1 ). 

These estimates were obtained assuming risk heterogen- 
eity governed by parameters a and y previously estimated 
from global tuberculosis data [8]. Additionally, we ana- 
lyzed the sensitivity of estimated proportions of individuals 
in each class and of R 0 estimation to variation in these pa- 
rameters (Additional file 3 and Table 4). To avoid very low 
(i.e., close to 0) and very high (i.e., close to 1) values of a 
and y, which greatly affect estimations, we consider values 
between 0.005 and 0.995. This sensitivity analysis results 
in particularly large variations in susceptible and latently 
infected classes (Table 4). Considering the Portuguese 
population size, the estimated values indicate approxi- 
mately 11,000 to 14,000 people with primary TB infection, 
600,000 to 1,600,000 with latent infection, 1,200 patients 
under treatment and 500 undiagnosed actively infected 
persons. Rq varies from 1.116 to 2.421 (/3 varies from 
72.358 to 156.98 yrs -1 ). 



Trends of TB prevalence with transmission coefficient 
and R 0 

Assuming that systems (1) and (3) are at equilibrium, we 
can calculate how the proportion of infected individuals 
(/ for system (1) and I\ + I* 2 for system (3)) vary with 
transmission coefficient /? (Additional file 2). Figure 3 
represents the proportion of infectious individuals in 
system (1) (full line) and system (3) (dashed-line) for 
different values of fi and R 0 . The curves were obtained 



by considering the parameter values in Table 1 with the 
proportion y and the susceptibility factor a of the low- 
risk group fixed to the values estimated in [8]. The 
dotted-lines indicate the values of ft and Ro for systems 
(1) and (3) corresponding to the estimated proportion of 
infectious individuals 4.261 x 10~ 5 . 

Discussion 

With a TB prevalence of over 20 cases per 100,000 indi- 
viduals, Portugal is one of the countries with the highest 
TB burden in the European Union (along with Bulgaria, 
Estonia, Latvia, Lithuania, Poland and Romania) [9]. 
Furthermore, there are worries that the current eco- 
nomic crisis may increase TB burden in this country 
due to a consequent increase of unemployment, wealth 
inequality and general deterioration of the living con- 
ditions [17,18]. For this reason, developing a mathem- 
atical model that suits the Portuguese scenario is 
important as to better understand the dynamics of this 
transmission system, identify key parameters in the 
control of TB and to inform policies by predicting the 
impact of interventions. 

Despite the clear need for a better understanding of TB 
in Portugal, this work is the first that considers the com- 
plexity of TB transmission to estimate the Portuguese 
TB burden. Its main limitations are: the calculation of 
TB detection time from retrospective self-diagnosed 
data; the simplifications of the model, particularly, not 
assuming age-structure nor emigration; and the as- 
sumption of a steady-state for TB transmission. As 
such, although the results provide a valuable indication 
of TB burden in each stage of TB infection, they should 
be interpreted with caution. 

In the present work, we had access to Portuguese 
clinical and socio-demographic data and calculated 
directly the standard detection time for TB as of about 



Table 4 Estimates of active disease prevalence and basic reproduction number for tuberculosis transmission models (1) 
and (3) 



System 3 


S 


P (10" 3 ) 


1 (10" 5 ) 


L 


T (10" 4 ) 


Ro 


1 


0.844 


1.059 


4.261 b 


0.155 


1.148 


1.116 


3 


0.930 


1.316 


4.261 b 


0.069 


1.148 


2.237 


3 with varying a and y 


pr: 0.844 


pr: 1 .059 


4.261 b 


pr: 0.059 


1.1 48 c 


pr: 1.116 




Pr: 0.940 


Pr: 1 .348 




Pr: 0.155 




Pr: 2.421 



Abbreviations: S Susceptible, P Primary infection, I Active Infection, L Latent Infection, R 0 Reproduction number, pr 0.5% percentile, Pr 99.5 
a - parameter values according to Table 1 unless stated otherwise. 
b - estimated directly from tuberculosis incidence. 
c - estimation is not a function of a or y. 



i percentile. 
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Basic reproduction number, R 
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Figure 3 Prevalence of active TB as a function of transmission 
parameters. Proportion of infectious individuals as a function of $ 
and fi 0 . Curves represent the endemic equilibria according to the 
homogeneous system (1) (full line) and the heterogeneous system 
(3) (dashed line), using parameter values specified in Table 1. 
Estimates of /3 and ft 0 for the same proportion of infectious 
individuals under both models are marked (dotted lines). 



2.81 months and the typical treatment length as 8.80 
months. These values are comfortably within the range 
recommended by WHO [1]. The default proportion in 
Portugal was calculated as 4%, a value that is somewhat 
lower than previously reported (6 to 30%) [19]. Using 
these calculated values along with parameters taken 
from the literature, we assessed our proposed model 
against the Portuguese data by comparing expected 
temporal dynamics with those recorded in the database 
(Figure 2). Agreement is impressive given the simple 
and unconstrained nature of the model, indicating its 
adequacy to the study of TB transmission in Portugal. 
For this reason, although more complex heterogeneous 
models have been proposed [20,21], we are confident 
on the suitableness of the considered model. 

The burden of TB in a country is difficult to measure 
due to the different stages of the disease and the impos- 
sibility to quantify directly from surveillance data the 
proportion of individuals in some of these stages. In the 
adopted model we considered four stages of TB infec- 
tion: primary infection; active disease not yet diagnosed; 
active disease undergoing treatment; and latent infec- 
tion. Surveillance data is limited to patients entering the 
treatment stage. However, by considering a specific 
mathematical model we were able to estimate the total 
burden of TB in Portugal as of about 1.6 million infected 
individuals (i.e., about 11,000 with primary TB infection, 
1,600,000 with latent infection, 1,200 under-treatment 



and about 500 undiagnosed active infections), correspond- 
ing to an R 0 of 1.116. This value for R 0 is very much in line 
with previous studies [22,23]. 

Nonetheless, we recognize that, as evidenced by Murphy 
et al. [24] and Gomes et al. [8], individual variation in 
susceptibility to infection may change dramatically the 
transmission dynamics of TB leading to significant dif- 
ferences in the interpretation of population measures. 
Murphy and co-authors consider genetic factors as a 
partial cause for heterogeneity, but, as previously noted, 
there is much uncertainty associated to these factors [25]. 
Another possible cause for heterogeneity is whether indi- 
viduals are vaccinated or not, however, given the inconsist- 
ent protection conferred by the Bacillus Calmette-Guerin 
vaccine, there is also considerable uncertainty related to 
the impact of vaccination [26]. Bacaer and co-authors [27] 
considered a heterogeneous model to study jointly TB and 
HIV epidemiology. This model divided the population in 
high and low risk of TB infection, depending if the persons 
were HIV positive or negative, and was developed to 
analyze data from South Africa, a country where about 
12% of the population is infected with HIV [28]. Dowdy 
et al. [29] considered a heterogeneous transmission model 
to study TB data from Rio de Janeiro. In this study the 
high and low risk division was done by considering geo- 
graphically distributed prevalence data, however, the ex- 
istence of clearly defined "hotspots" of TB cases in this 
setting is mostly due to slumbers and other well-defined 
regions with a high impact of risk factors and social de- 
terminants. We argue, then, that heterogeneity in TB 
infection risk in Portugal is likely to result from a differen- 
tial impact of biological and social determinants to differ- 
ent sections of the population [3,17]. An often neglected 
caveat of introducing heterogeneity in TB transmission 
models is that the factors that drive heterogeneity may 
break assumptions of random mixing (e.g. being in prison 
or nursing-home, working on health-care, living in slum- 
bers). In Portugal, however, the factors that have been im- 
plied with heterogeneity in TB risk are not markedly 
associated with closed-groups (e.g. HIV and drug-use [9], 
unemployment, wealth inequality and living conditions 
[17]). Given the difficulty to categorize the Portuguese 
population as having a high or low risk of infection so 
that the world-wide estimations from [8] could be dir- 
ectly verified and refined, we calculated the TB burden 
from a spectrum of profiles of heterogeneous suscepti- 
bility. Hence, we estimated a range of possible TB bur- 
den of between 0.6 and 1.6 million infected individuals 
(i.e., 11,000 to 14,000 with a primary TB infection, 
600,000 to 1,600,000 latent-infected, and the same 1,200 
under-treatment and 500 undiagnosed actively-infected 
as verified in the homogeneous scenario), correspond- 
ing to R 0 from 1.116 to more than 2 All. The wide 
range of values for TB burden is mostly attributed to 
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the wide range of estimates of the proportion of latently- 
infected individuals. Considering an increasing trend on 
the reactivation rate of latent infections [13], it is of great 
importance to characterize the heterogeneity of a popula- 
tion and, thus, to estimate with accuracy the proportion of 
latently-infected individuals in a population. 

In the comparison between homogeneous and hetero- 
geneous models, although we considered the same pro- 
portion of infectious individuals (i.e. 4.261 x 10~ 5 ), we 
estimated a possible range of values for R 0 and /? in the 
heterogeneous case that can be considerable distinct 
from the ones estimated in the homogeneous case. This 
interesting dynamic is better explored in Figure 3, 
where we calculated the proportion of infected individ- 
uals for various values of /? and R 0 while considering 
systems 1 and 3 (with fixed values of a and y). These re- 
sults show that TB estimates of transmission coefficient 
/? (and of R 0 ) may vary greatly depending on the model 
considered. Examining the effect of heterogeneity in R 0 
is especially important since, as noted before, its use in 
TB may be problematic if all the dynamics of the trans- 
mission system are not considered [30,31]. 

Conclusions 

We proposed a mathematical model for TB transmis- 
sion, which fits well the dynamics described by records 
of the Portuguese National Directorate of Health, as a 
tool for intervention impact studies. Under homogen- 
eity assumptions, we estimated an upper bound for the 
total burden of TB infection in Portugal of 1.6 million 
individuals (more than 15% of the Portuguese population). 
Accounting for heterogeneities in susceptibility and ex- 
posure to infection led to lower estimates of infection 
prevalence. 
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